---
abstractspacing: double
fontsize: 12pt
margin: 2cm
urlcolor: darkblue
linkcolor: Mahogany
citecolor: Mahogany
spacing: single
bibliography: references.bib
biblio-style: apalike
output:
  pdf_document:
    citation_package: natbib
    fig_caption: no
    number_sections: no
    keep_tex: no
    toc: no
    toc_depth: 3
    template: article-template.latex
---

```{r setup, include=FALSE}
knitr::opts_chunk$set(echo = TRUE, warning=FALSE, message=FALSE,cache=TRUE)
knitr::opts_chunk$set(fig.width=7, fig.height=6, out.width = '70%', fig.align = "center") 
rm(list=ls())
library(remotes)
library(kableExtra)
library(haven)
library(tidyverse)
library(ivmodel)
library(doParallel)
library(foreach)
library(estimatr)
require(AER)
library(lfe)
library(glue)
#path <- "/Users/ziwenzu/Dropbox/research/IV/IV Sensitivity/LLXZ_rep"
path <- "~/Dropbox/ProjectZ/IV Sensitivity/LLXZ_rep"
setwd(path)
knitr::opts_knit$set(root.dir = path)
#install_github("apoorvalal/ivDiag")
library(ivDiag)
# number of cores
cores <- 15
```


## @trounstine2016

| Replication Summary | |
|----------------|-----------------------------|
| Unit of analysis | city\*year |
| Treatment | racial segregation |
| Instrument | the number of waterways in a city; logged population |
| Outcome | direct general expenditures |
| Model | Table5(1)|

```{r ajps_Trounstine_2016}
df<-readRDS("./rawdata/ajps_Trounstine_2016.rds")
D <-"H_citytract_NHW_i"
Y <- "dgepercap_cpi"
Z <- c("total_rivs_all", "logpop")
controls <- c("dgepercap_cpilag","diversityinterp","pctblkpopinterp",
              "pctasianpopinterp","pctlatinopopinterp","medincinterp",
              "pctlocalgovworker_100","pctrentersinterp","pctover65",
              "pctcollegegradinterp","northeast","south","midwest", 
              "y5", "y6", "y7", "y8", "y9")
cl <- NULL
FE <- NULL
weights<-NULL
(g<-ivDiag(data=df, Y=Y, D=D, Z=Z, controls=controls, FE =FE, 
  cl =cl,weights=weights, cores = cores))
```

```{r, cache = FALSE}
plot_coef(g)
```


```{r ajps_Trounstine_2016_sav, echo = FALSE}
save(g, file="./estimate/Trounstine2016.RData")
```


